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Abstract 

Composite pulses have found widespread use in both conventional Nuclear Magnetic Resonance experiments and in 
experimental quantum information processing to reduce the effects of systematic errors. Here we describe several 
families of time symmetric and antisymmetric fully compensating composite pulses, inspired by the previous F„, G„ and 
BBl families family developed by Wimperis. We describe families of composite 180° pulses (not gates) which exhibit 
"unprecedented tolerance of pulse strength errors without unreasonable sensitivity to off-resonance errors, and related 
families with more exotic tailored responses. Next we address the problem of extending these methods to other rotation 
"angles, and discuss numerical results for 90° pulses. Finally we demonstrate the performance of some 90° and 180° 
.pulses in NMR experiments. 



Q-il. Introduction. 
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NMR experiments suffer from a range of errors which 
can be traced back to imperfections in the radiofrequency 
P^HiControl fields used to manipulate spin systems. In par- 
—iticular the Bi field is not infinitely strong, leading to off- 
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resonance errors, and is not uniform across the sample, 
■leading to pulse strength errors. One particularly success- 
l/~, ful approach for tackling such errors is the use of com- 
(^ posite pulses [ij, l2| , in which a single rotation is replaced 
C^ ^jy a series of rotations chosen such that in the absence 
CO of errors the combined propagator implements the desired 
(T^ rotation, while in the presence of small errors the effects 
^^ of the errors in individual rotations largely cancel one an- 
C^ pther. Here we principally consider the problem of tack- 
. . ding pulse strength errors without introducing additional 
^ Sensitivity to off-resonance errors. 

Corresponding errors arise in experimental implemen- 
Itations of quantum information processing (QIP) [3[ |j| 
Ci where they arc known as systematic errors [5|, to distin- 
guish them from random errors arising from decoherence, 
and where there is considerable interest in performing ex- 
tremely accurate unitary transformations on quantum sys- 
tems in the presence of realistic experimental errors. Again 
composite pulses provide a potential solution, but there are 
many significant differences between the design of compos- 
ite pulses for quantum computing and conventional NMR. 
Firstly, many traditional composite pulses, are designed 
to act correctly only on particular initial states, and so are 
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not suitable for quantum computing, where fully compen- 
sating (Class A) composite pulses '2i] have to be used [5|- 
llSJ I . (While fully compensating pulses can be used for con- 
ventional NMR experiments, where they have the advan- 
tage that they can be immediately substituted for any sim- 
ple pulse without the need for detailed analysis '^, their 
use is normally excessive and better results can normally 
be obtained using pulses tailored to the relevant problem.) 
Similarly it is necessary to use genuinely fully compensat- 
ing pulses, rather than using phase cycling or gradient co- 
herence pathway selection to suppress imperfections 14| . 



Secondly, although composite pulses to suppress off- 
resonance errors, originally developed by Tvcko |15l| , have 
been used for QIP in NMR P, [^ and SQUID (Tjexperi- 
ments, tackling off-resonance errors is (with the exception 
of dynamic decoupling, discussed below) rarely of much 
interest in quantum computing, as the aim is to control 
a previously well characterised spin system, rather than 
to investigate an unknown system. It is, however, gener- 
ally desirable to ensure that the sensitivity to off-resonance 
errors is not greatly increased, and in particular is not sig- 
nificantly worse than that of a simple pulse. 

Thirdly, it is useful to distinguish between composite 
pulses designed for single use, and pulses designed for use 
in a pulse train, such as decoupling sequences, where it 
may be desirable to combine relatively simple compos- 
ite pulses with phase cycles and supercycles IT], [lE 



In 



quantum computing most work has concentrated on sin- 

30, 



gle pulses, but in the field of dynamic decoupling 
which is likely to play an important role in the design of 
quantum memories, the use of techniques adapted from 
conventional decoupling sequences has proved fruitful 21 , 
|22| . Here we will only consider the case of single pulses. 
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Finally, quantum computing generally requires very ac- 
curate control in the presence of small or moderate errors, 
rather than moderately accurate control with very large er- 
rors. (These two effects sometimes go together, but some- 
times do not). For this reason it is common to express the 
fidelity of the pulse as a Taylor series as described below, 
and to seek to suppress as many error terms as possible. 
In particular there has been significant interest in the pos- 
sibility of develo p in g families of arbitrarily accurate com- 
posite pulses 23|, |2J1 , although actually designing pulses 
by such methods is normally a complex problem [25|. 

Here we will show that it is straightforward to find ro- 
bust 180° pulses with arbitrary tolerance of pulse strength 
errors, but extending these ideas to pulses with other rota- 
tion angles is more challenging and only partially solved. 
As with other composite pulses developed in NMR these 
can be extended to apply to spin-spin couplings [26[ and 
can also be applied in a wide range of other quantum com- 
puting experiments |27H3C| . We will restrict ourselves to 
9o rotations about the x-axis of the Bloch sphere; rota- 
tions about other axes can be trivially derived from these 
by offsetting all pulse phases by the desired phase angle. 

2. Pulse fidelity 

It is convenient to characterize simple and composite 
pulses by their propagator fidelity [9| 



The BBl pulse sequence, originally developed by Wim- 



peris [33|, is effective in such cases; it replaces a ^o pulse 



T= |tr(y[/^)|/2 



(1) 



where V is the propagator of the pulse in the presence of 
errors and U is the propagator for the ideal pulse (tak- 
ing the absolute value is necessary in general to neglect 
the effects of global phases, and dividing by two simply 
normalises the fidelity). Equivalently pulses can be cate- 
gorized by their infidelity, defined hy X — 1 — J^. 

Much work to date has concentrated on pulse strength 
errors, which arise from errors in the strength of the Bi 
field. These are closely related to pulse length errors, and 
this name is often used instead, but a careful distinction 
between the two errors can and should be made 3l|. In 
particular the two cases are not quite the same in the pres- 
ence of simultaneous off-resonance effects, as it is impor- 
tant to distinguish between measuring the off-resonance 
field as a fraction of the nominal Bi field and as a fraction 
of the actual Bi field S^- In reality both types of error can 
be present: while naturally occurring errors usually arise 
from inhomogeneity in the Bi field strength, many exper- 
imental demonstrations of composite pulses introduce ad- 
ditional deliberate errors by varying the pulse length [31| . 

In the absence of off-resonance errors, pulse strength 
errors replace the rotation angle 6 of each pulse by 0' = 
(1 + e) 9, where e is the fractional error. The pulse fidelity 
can then be conveniently expanded as a Taylor series in e; 
for a simple pulse this leads to 



by the pulse sequence 



TTp 27r3/3 TT^Oo /3 = ± arccos(— ^/47r) 



(3) 



written with time running from left to right. (The choice 
of sign is arbitrary but must be made consistently.) The 
initial three pulse sequence performs no overall rotation in 
the absence of errors, but in the presence of pulse length 
errors generates a pure error term which largely cancels 
the error in the final 6o pulse, suppressing both second 
and fourth order infidelities. While Wimperis originally 
applied the correction sequence at the start of the com- 
posite pulse, as shown above, it can instead be placed at 
the end of the sequence, or indeed in the center of the 
sequence by splitting the 0o pulse into two equal halves 
[g]. We wiU initially only consider the case of 180° pulses; 
these are used in NMR QIP to implement not gates and 
in more conventional NMR to implement spin echoes and 
their generalisation dynamic decoupling 
delity of a BBl 180° pulse is 
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The fi- 



J'bbi = l-e^ X 5 7rVl024+O(e**). 



(4) 



T = I cos(e6'/2)| = 1 - e^e^/8 + O(e^). 



(2) 



Note that as usual a propagator infidelity of order 2n cor- 
responds to an error term in the underlying propagator of 
order n [9|, and so the sixth-order sequence BBl removes 
the first and second order errors from the propagator. 

BBl has found widespread use in quantum comput- 
ing experiments, but there is interest in designing Class A 
pulses with even greater error suppression. (While the ap- 
proaches of arbitrary precision composite pulses have been 
used to design candidate pulses [23|-|25| , those explored to 
date have not proved useful as their theoretically superior 
performance is not confirmed in experiments [8J.) The 
most naive approach, constructing a higher order com- 
posite pulse by iteratively replacing each component pulse 
with a composite pulse, is rarely successful. This is not 
surprising, as in general there is no reason to expect that 
the residual errors in a composite pulse should have the 
same form as the errors in a simple pulse, and so it is un- 
likely that the very same approach will work at different 
orders of iteration. 

Instead methods for generating higher order pulses by 
recursive or iterative procedures (e.g., [32l l3J, |35|) nor- 
mally work directly on the error propagator 36]. An un- 
usual partial exception is the S„ family of composite in- 
version pulses [33]; although this family was designed by 
considering the error propagator, the resulting sequences 
are in fact simply naive iterative expansions. However the 
S„ family are not good Class A pulses, and so are not gen- 
erally suitable for quantum computing. More generally 
Tycko et al. have considered the use of naive iterative ex- 
pansions to design inversion and excitation pulses [37[ , and 
have discussed the conditions under which this approach 
is successful. 

As we will discuss below, various families of Class A 
composite pulses described by Wimperis, can be described 



by a simple iterative process, and thus provide a simple 
route to arbitrary precision 180° degree pulses for use in 
quantum computing. Such pulses are also moderately ro- 
bust at very large pulse strength errors, and so could in 
principle find applications in experiments using very in- 
homogeneous Bi fields, such as phase-modulated rotating- 
frame imaging jSStj or single-sided NMR [39|], as well as 
effectively removing the need to calibrate RF pulse widths 
in walkup NMR experiments. However in conventional 
NMR experiments it is likely that tailored approaches will 
be more effective. While many suitable pulses have been 
developed using conventional approaches, we also draw 
attention to the robust inversion sequences developed by 
Torosov [40|, while in walkup NMR shaped pulses devel- 
oped u sing optimal control theory 4l| , such as the Fanta4 
pulses [43 may also be of interest. 

3. Antisymmetric composite pulses 

The Wimperis error correcting sequence is time sym- 
metric, and the BBl sequence as a whole will similarly be 
symmetric if the correction sequence is placed in the mid- 
dle of the main pulse. As we will see later such symmetric 
pulses have certain advantages, especially in the presence 
of off-resonance errors, and many composite pulses have 
this property. Here we instead consider time antisymmet- 
ric composite pulses 180° pulses, which can have partic- 
ular advantages in some conventional NMR experiments 



14 l43l |. As we shall see, it is possible to combine many 
of the best features of the two approaches, by first de- 
signing antisymmetric composite pulses and subsequently 
converting these into corresponding symmetric pulses, or 
vice-versa. 

The design of antisymmetric composite tt pulses rests 
on the observation that, in the absence of errors, any time 
antisymmetric sequence of tt pulses such as 



■ 7r_0i TTo TT^j 



■^0„ 



(5) 



is equivalent to a simple ttq pulse [43|. Thus the phases of 
the outer pulses, {0i, 02, • ■ • 4>n\, can be adjusted at will, 
with the possibility of creating error tolerance. For the 
case of a five pulse sequence, with two controllable phase 
angles, the system can be solved analytically, and the 
second order infidelity term can be removed by choosing 
01 = ±0 and 02 = ±30 with = arccos(-l/4) « 104.5°. 
Once again the choice of sign is arbitrary, but must be 
made consistently, and in this case negating the phases 
corresponds to replacing the composite pulse with its time 
reversed form. In passing we note that such antisymmet- 
ric pulses have the useful property that the rotation axis 
of the overall propagator always lies in the xz plane [37[ , 
a fact which has considerable importance in some conven- 
tional NMR applications [145, but which is not directly 
relevant in QIP. 

Intriguingly these phases are identical to the BBl phase 
angles for a tt pulse, and just like BBl this pulse also 



suppresses the fourth order infidelity; indeed the fidelity 
of the antisymmetric pulse is identical to that of the BBl 
TT pulse. This is not coincidental, but rather can be traced 
to two simple facts. Firstly the fidelity definition, Eq. [U 
is invariant under cyclic reordering, and secondly we can 
use the identity 

TTo 6*0 = fl.^TTo (6) 

for a ttq pulse and any other pulse. Thus 

•^ = tr (tt^^, tt^^^ tt;, tt^^ tt^^ ttq) /2 (7) 

= tr (vr^ vr^^ tt^^ tto tt^^^ tt^^i ) /2 (8) 

= tr{iT'^n'^^n'^^T:'^^T:'^^T:o)/2 (9) 

where tt' — 7r(l+e) as before. Any antisymmetric sequence 
of five TT pulses can be converted into the same form as a 
BBl pulse, and so the optimal error correcting sequences 
will occur at the same phase angles. The ambiguity in 
the sign of the arccos function in Eq. [3] is now seen to 
correspond to the equivalence of an antisymmetric pulse 
sequence and its time-reversed form. 

This antisymmetric composite pulse was previously de- 
scribed by Wimperis [ij, |43|, who called it an Fi pulse. 
It was originally derived by using a toggling frame de- 
scription of an antisymmetric composite pulse to isolate 
the error term, and then removing this term to first or- 
der J43| . Fi is, however, simply the first member of the 
F„ family of pulses with rapidly growing error tolerance. 
Wimperis explicitly described F2, the second member of 
the family, and gave an outline prescription for designing 
further members. 

Wimperis designed the family using an iterative expan- 
sion of the error propagator, but it is instead possible to 
proceed by a naive iterative expansion, effectively replac- 
ing each TT pulse in an Fi sequence by an Fi sequence. 
It is, however, necessary to replace alternate pulses by 
"standard" Fi sequences and their time-reversed equiva- 
lents. This naive iterative approach to the design of Class 
A composite pulses is not usually successful, but its appli- 
cation to time antisymmetric sequences of tt pulses, and 
the need to alternate standard and time-reversed variants, 
can be understood h^ considering the composite pulse in 
the toggling frame [43| . 



By contrast with other approaches [23l - l25j the F„ fam- 
ily provides a simple and direct method for explicitly con- 
structing Class A composite pulses with arbitrary preci- 
sion: the whole family can be easily described by listing 
the phases of successive tt pulses, with the iterative form 

0n+l = {-30 + 0„, -0 - 0„, 0„, - 0„, 30 -|~ 4>n} (10) 

with 00 = and — arccos(— 1/4). Thus 

0i = {-3, -1, 0, 1, 3}x0 (11) 

reproducing Fi, and 

02 = {-6, -4, -3, -2, 0, 2, 0, -1, -2, -4, -3, 

- 1, 0, 1, 3, 4, 2, 1, 0, -2, 0, 2, 3, 4, 6} x (12) 



pulses 






1 


1 


5 


2 


25 


3 


125 


4 


625 



order 
-71 



,18 



^54 
,162 



coefficient 

"^5^72^ 

5 X ^6/2'° 
55 X 7riV23i 
513 

540 



numerical value 



5^^ X 7r''V29 



X TT 



162 /o283 



72^ 



1.234 
4.694 
258.6 

4.324 X 10^ 
2.021 X 10^3 



Table 1: Summary of the performance of F„ pulses for moderate 
values of n. F„ pulses are only perfect at e = but near this point 
show infidelities of the orders shown. 



describes the 25 pulses making up F2. In this notation 
the direct implementation of a not gate as a single pulse 
can be considered as the family member Fq. Note that 
the alternating pattern of plus and minus signs in Eq. [TU] 
corresponds to the alternate use of standard and negative 
phase (time-reversed) pulses; this is necessary because the 
phase of the error term in the toggling frame is alternately 
positive and negative [43[ , and the phase steps of the com- 
ponent pulses must be alternately negated to follow this 
pattern. 

The F2 composite pulse suppresses all infidelity terms 
below eighteenth order, with 



J'f, 



1 



e^^ X 625 TT 



18 



/2147483648-f 0(e 



20 \ 



(13) 



and the performance continues to improve rapidly for higher 
members of the family, as listed in Table [TJ (These results 
were obtained by direct calculation with pulse sequences 
of the form of Eq. (TU]) In general F„ is a sequence of 5" 
pulses with a fidelity given empirically by 

J-p^^ = 1 _ £29 X St'?-!)/^ X TT^? X 2(l-7«)/2 +0(e29+2) (14) 

with q — 3" (this expression has not been checked be- 
yond F5, which suppresses all infidelity terms below 486*'' 
order). For the higher members of the family the first non- 
zero coefficient is very large, but has little effect at mod- 
erate e as it is combined with an extremely high power of 
e. Note that when evaluating the properties of sequences 
beyond F2 it is essential to use either fully analytic meth- 
ods or very high precision numerical calculations in order 
to capture fully the complex pattern of cancelations that 
leads to very high order error tolerances. 

The performance of some of the lower members of the 
F„ family is shown in Fig.[l] clearly showing the very broad 
error tolerance achievable. Theoretically these pulses al- 
low astonishingly precise rotations to be performed, but 
in practice such extremely precise rotations are unlikely to 
be cither required or achievable, and for this reason the in- 
fidelity plots shown here concentrate on infidelities around 
10~^. These practical limits arise because no experimen- 
tal implementation will ever be completely perfect in every 
regard. In particular, it is implausible that the required 
phase shifts can be produced precisely, and it is likely that 
relaxation during pulses will become a problem with the 
very long pulses required for high members of the family. 

A similar approach can be used to describe the G„ 
family of composite pulses [43|. Rather than seeking to 




Figure 1: Fidelity T and infidelity X as a function of pulse strength 
error e for the F„ family of pulses from Fq (plotted in red) to F3 
(plotted in black). 



create a composite pulse with high fidelity over a central 
region, these pulses are designed to give perfect fidelity 
at certain particular error values, in the hope that this 
will lead to moderate error tolerance over a wide range of 
intermediate values. (Error tolerance of this kind is not 
likely to be particularly useful for quantum information 
processing, but may find applications in more conventional 
NMR.) 

The first non-trivial member, Gi, is obtained by form- 
ing an antisymmetric sequence of tt pulses as before and 
then choosing the phases to maximise the pulse fidelity at 
e = ±0.5. As before the problem can be solved analyti- 
cally, giving the sequences of phases 



7i = {l, -2, 0, 2, -1}X7 



(15) 



with 7 = 7r/4, a sequence which Wimperis called Gi. The 
antisymmetric structure of the pulse guarantees that it will 
also be perfect at e = 0, and the infidelity is quadratic in e 
around these fixed points, with a shallow minimum in the 
fidelity around e w 0.28. Converting this sequence into a 
time symmetric form gives the Wimperis composite pulse 

BB2 [3aj- 

There is nothing special about the choice of e = ±0.5 
for the points of ideal fidelity, and close relatives of Gi 
can be obtained for any value up to 0.8. However beyond 
0.5 the depth of the fidelity minimum increases rapidly, 
removing the desired broad error tolerance, while for small 
values the composite pulse sequence becomes very similar 
to Fi. We will, therefore, stick to the original choice of 
0.5. Once again this composite pulse can be expanded 
iteratively to get higher members of the G„ family. The 
rule for pulse phases is entirely analogous to that for F„, 
and is 

7n+l = {7 + In, -27 - 7n, 7n, 27 - 7„ , -7 + 7„} (16) 



pulses 






1 


0.000 








1 


5 


0.000 


±0.500 






2 


25 


0.000 


±0.500 


±0.786 




3 


125 


0.000 


±0.500 


±0.786 


±0.911 


4 


625 


0.000 


±0.500 


±0.786 


±0.911 



±0.963 



Table 2: Gn pulses exhibit perfect fidelity at 2n + 1 points with the 
approximate values shown. 





Figure 2: Fidelity T and infidelity I as a function of pulse strength 
error e for the Gn family of pulses from Go (red) to G3 (black). Note 
that only the middle three infidelity minima occur at rational values 
of e, and the outer minima are not well digitised. 



with 7o = and 7 = 7r/4. G2 has perfect fidehty at e = 
and e = 0.5, just Hke Gi, but also has perfect fidelity at 
two additional points at e « ±0.786. In general G„ has all 
the same perfect points as G„_i, and a pair of additional 
points further out, as shown in Table [5] and Fig. [5] Beyond 
Gi there is no simple formula for the location of these 
points, which must be found numerically. 

4. Narrowband and passband pulses 

In addition to describing two symmetric families of 
broadband sequences, BBl which corresponds to Fi and 
BB2 which corresponds to Gi, Wimperis also described 
33] two families of narrowband sequences (which only give 
effective pulses when e is close to zero) and two families 
of passband sequences (which act as error correcting com- 
posite pulses for moderate values of e, but act as mod- 
erately robust IDENTITY Operators for values of e near 
±1). As usual we will consider the members correspond- 
ing to a nominal ttq rotation. These families are all time- 
symmetric, but as before we can convert them into time- 
antisymmetric sequences of tt pulses. Having done this we 
can use the iterative approach to form higher members of 



Figure 3: Fidelity J^ as a function of pulse strength error e for the Nn 
family of pulses from Np (plotted in red) to N3 (plotted in black). 



the same family, and then (if desired) reverse the transfor- 
mation to create time-symmetric sequences. 

We will only discuss the narrowband family derived 
from NBl, which we call N„, and the passband family 
derived from PBl, which we call P„. Applications of pass- 
band pulses in solid state NMR have recently been de- 
scribed [4J], and we will demonstrate an application of 
narrowband pulses below. 

An NBl ttq pulse takes the form 



"Kp 27r_i, TT^ TTo V — arccos(— 1/4) 



(17) 



and can be rearranged to get the antisymmetric pulse Ni, 
defined by the phases 



v\ 



{1, -1, 0, 1, -1} X V. 



(18) 



(A related pulse was previously discussed by Tycko et al. 
[37[ , but they only considered its properties as an inversion 
pulse.) As usual the higher members are defined by 



'^n+l 



{v+Vn, 



-V - l/„, Z/„., V 



-P+Vn}. (19) 



All these pulses have second order infidelity around e = 0, 
with the coefficient of this error term rising with increasing 
n 

J"n„ = 1 - e^ X ttVS X (15/4)" + 0{e^), (20) 

creating a pulse which is only effective over increasingly 
narrow regions, as shown in Fig. |31 
A PBl TTp pulse takes the form 



27r^ 47r_^ 27r^ ttq ip = arccos(— 1/8) 
with a fidelity compared with an ideal n pulse of 

J'pBi = 1 - e^ X 63 7rVl024 + 0{e^) 
and a fidelity compared with an identity gate of 
J- = 1 - (1 + e)'^ X 63 7rV512 + 0{{l + ef) 



(21) 



(22) 



(23) 



where the expansion of the fidelity has been taken around 
e = — 1. PBl has slightly worse suppression of small errors 
than BBl, but gives very little excitation in regions of 
extreme pulse strength error. 




0.5 e 1.0 




Figure 4: Fidelity as a function of pulse strength error e for the P„ 
family of eomposite pulses from Pq (red) to P3 (black). 



PBl can be rearranged to give the Pi sequence 



7r_^ TT-^ n^ TT^ TTq TT-,p TT-^ TT^ TT^ 



(24) 



which has exactly the same fidehty but is now described 
as a time-antisymmetric sequence of 9 tt pulses. (This 
pulse was previous discussed by Odedra and Wimperis [4J| 
under the name APBi.) This can be expanded iteratively 
to obtain P2 and so on. The P2 sequence of 81 tt pulses 
has a fidelity 



-^P. 



1 



,18 



x3«7V^V2^^ + 0(e^") 



(25) 



with the expected eighteenth order behaviour. However 
the fidelity to an identity gate around e = — 1 remains 
fourth order, with the same coefficient. The behaviour of 
higher members of the family is shown in Fig. 21 demon- 
strating increasingly robust passband behaviour. 

5. Combining pulses 

The results for F and G pulses discussed above are 
largely implicit in the earlier work of Wimperis [43[, al- 
though he did not give an explicit form for F„ and G„, 
or derive the remarkable error tolerances shown here, but 
the results for N„ are new, as are those for P„ for the 
sequences beyond Pi. Our approach is not, however, con- 
fined to these four families. 

All these composite pulses are built up iteratively, by 
replacing each n pulse in a sequence by an error corrected 
TT pulse, but it is not necessary to use the same error- 
corrected pulse at each stage of the iteration, and the 
phase patterns described above can be applied more gen- 
erally. (Similar ideas have been explored in the design of 
composite inversion and excitation pulses |37l.l45|.) In par- 
ticular, consider the GF pulse obtained by applying the G 



Figure 5: Fidelity J-" and infidelity X as a function of pulse strength 
error e for F (red), GF (green), FGF (blue) and FFGF (black) com- 
posite pulses. 



pattern of phases to Fi 

i>GF = {7 + (/)i, -27 - (/)i, 01, 27 - (/)i, -7 + 0i} (26) 

to give a composite pulse made up from 25 n pulses. In 
common with Fi this pulse has sixth order infidelity around 
e = 0, and in common with Gi this pulse has perfect fi- 
delity at two additional values of e, although in this case 
they are moved out beyond 0.5, with quadratic error terms 
around these additional points. This results in good error 
tolerance over the central region, and moderate error tol- 
erance over a much wider region as shown in Fig. [5l These 
new perfect points occur at e « ±0.720 where the underly- 
ing Fi pulse has a fidelity of l/\/2, the same as the fidelity 
of a naive pulse at e = 0.5. 

In the same way the FG pulse can be obtained by ap- 
plying the F pattern of phases to Gi 

i'FG = {-30 + 711 -</>- 71, 71, 0-71, 30 -I- 71} (27) 

to get another 25 pulse sequence. In common with Gi this 
pulse exhibits perfect fidelity at e = ±0.5, and in common 
with Fi it has sixth order infidelity around e = 0. However 
this pulse is also sixth order around e — ±0.5, giving good 
error tolerance over a wide region. 

The process can be extended in the obvious way: for 
example the 125 pulse sequence FFG (or F2G), obtained 
by applying the F phase pattern to FG, has perfect fidelity 
at e = and e — ±0.5, and then shows eighteenth order 
infidelity around these three points. Similarly the sequence 
FGF has eighteenth order infidelity around e = and sixth 
order infidelity around e w 0.720 as shown in Fig. [51 Any 
combination of F and G iterations can be applied: each 
application of G creates an additional pair of points at 
which the composite pulse has perfect fidelity, and each 




Figure 6: Fidelity J-' as a function of pulse strength error e for a 
G pulse (red), NG (green), N2G (blue), and N3G (black); the N3G 
pulse is only effective near e = and e = ±0.5. 



application of F triples the order of the infidelity around 
any existing perfect points. 

While the most obvious applications of this approach 
lie in the creation of broadband pulses, such as FGF, it 
is also possible to consider more exotic combinations. For 
example, applying the N pattern of phases to a G pulse 
gives a family of composite pulse which in effect select for 
three particular control field strengths. All such pulses 
show perfect fidelity at e = and e — ±0.5, but as the N 
pattern of phases is repeatedly applied the fidelity away 
from these points is reduced, as shown in Fig. |6l 

6. Off resonance errors 

This ability to sculpt the tolerance of pulse strength 
errors is remarkable, but it is important to check that it 
is not bought at the cost of increased sensitivity to other 
errors, in particular off-resonance errors. Note that we 
do not consider here the design of pulses which are si- 
multaneously resistant to off-resonance and pulse-length 
errors; discussions of this problem can be found elsewhere 

In the presence of simultaneous errors the propagator 
for a nominal 9^ pulse is 



V = exp(— i6'[(l + e){ax cos(j> + 1 



+ M]/2) (28) 



where e is the pulse strength error as before, and / is the 
off-resonance fraction defined here as the ratio between the 
frequency offset of the transition from the pulse and the 
nominal rotation frequency of the pulse in the absence of 
any errors. (Minor differences between the results shown 
here and those in previous work [6| reflect the fact that 
this earlier work considered pulse length errors rather than 
pulse strength errors.) In the absence of pulse strength 
errors the pulse fidelity is given by 



F ^l 



/2 X sin^ 



{e/2)/2 + 0{f' 



(29) 



so the pulse infidelity is largest when = tt, for which 
I K, /^/2, and smallest when 6 = 2mT, where the first 
order error term completely cancels out. 

In the absence of pulse length errors the BBl compos- 
ite pulse has the same sensitivity to small off-resonance er- 
rors as a single pulse; this favorable behavior can be traced 



to the fact that the time symmetric correction sequence is 
simply a 27r pulse inside another 27r pulse, and such nested 
sequences of 27r pulses have no first order error terms, so 
that the dominant source of infidelity is the main Oq pulse. 
This is not true for F„ and G„ pulses: direct calculations 
indicate that each iteration of the F phase sequence in- 
creases the size of the quadratic infidelity term by a factor 
of 16, while each iteration of G increases the infidelity by 
9 -I- 2-\/2 « 11.8. Thus for an FG pulse, for example, the 
off-resonance infidelity is X « 189.2/^. 

This greatly increased sensitivity to off-resonance er- 
rors might appear a critical flaw, as it would render the 
higher members of the F and G families completely unus- 
able in practice, but fortunately it can be easily overcome. 
As noted above the BBl composite tt pulse is very closely 
related to Fi, and can be obtained from it by moving the 
front part of the pulse sequence to the back with all the 
phases negated 



7r_30 7r_0 ttq tt^ n^^ 



7''o7r0 7r30 7r3^7r0. 



(30) 



As implemented above all the pulses before the main ttq 
pulse have been moved, but it is also possible to move 
this pulse as well, or to move half of the pulse to the 
back, resulting in a fully time-symmetric sequence, which 
has the advantage that the fidelity is then a purely even 
function of /. This transformation can be applied to any 
time antisymmetric sequence of n pulses, such as the F and 
G families, producing a nested sequence of 27r pulses: the 
extreme tolerance of pulse strength errors is retained, but 
now with little cost in sensitivity to off-resonance errors. 

The calculations above only apply in the absence of 
pulse strength errors, and it is necessary to consider the 
effects of simultaneous errors which lead to cross terms in 
the fidelity expression. These are best investigated numer- 
ically, as shown in Fig. [7] which plots the fidelity of various 
simple and composite pulses over a wide range of errors. 
While the simple FG sequence shows much greater toler- 
ance of pulse strength errors than BBl, this is achieved at 
the cost of greatly increased sensitivity to off-resonance er- 
rors. In contrast the time symmetric FG sequence achieves 
the same pulse strength tolerance as the antisymmetric FG 
sequence with only a modest increase in sensitivity to off- 
resonance errors compared with a simple pulse. 

7. Other pulse angles 

With the exception of the BBl pulse, Eq. [3l all the re- 
sults described above only apply to ttq pulses (not gates), 
and it would obviously be desirable to extend these to a 
wider range of pulse angles (the extension to tt pulses with 
other phase angles is trivially achieved by offsetting all the 
phases in a ttq composite pulse by the desired phase angle) . 
In particular it would be helpful to extend these methods 
to 7r/2 pulses, both for applications in conventional NMR 
and in QIP. Note that in QIP the set of 7r/2 pulses with 
arbitrary phases is universal for single qubit gates, and in 
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Figure 7; Fidelity as a function of off resonance error / and pulse strength error e for NOT gates implemented using simple pulses and BBl, 
FG, and symmetric BFG composite pulses; contours are drawn at fidelities of 0.9, 0.99 and 0.999. Note that the subsidiary maxima in the 
FG and symmetric FG plots are not normally useful, and the key indicator is the size of the principal maximum in the centre of the plot. 



NMR systems a robust 7r/2 rotation can be extended by 
analogy to produce a controUed-NOT gate [26|, and thus 
robust universal quantum logic. 

This process might seem simple, as the BBl composite 
pulse, Eq.[31 works for any rotation angle, and the BBl ttq 
pulse can be converted into an Fi pulse. Such optimism is, 
however, unjustified, as the identity in Eq. [6]only applies 
for TT pulses, and so the transformation can not be applied 
more generally. Furthermore the key property of antisym- 
metric pulses of the form of Eq. [5] only applies when the 
pulses are tt pulses. For example, the sequence 



bl f (701 f 02 



^^ 



(31) 



is only equivalent to ^o if </'2 = </'i + tt, and in general it 
is necessary to introduce similar relationships between the 
pulse phase angles. While we have been able to develop 
some low order antisymmetric 7r/2 pulses, for applications 
in QIP these sequences are inferior to the time symmetric 
sequences discussed below, and are not considered further 
here. 

We begin by rewriting the correction sequence at the 
start of the BBl pulse, Eq.[21 in the general time-symmetric 
form 

TT^i 7r02 7702 7^01 (32) 

which we refer to as a Wi sequence. Clearly such a pulse 
can be summarised by listing the two phase angles, (pi 
and ^2, which can be adjusted to optimise the fidelity of 
the combined pulse by removing terms in the Taylor series 
expansion. In particular it is possible to eliminate the 
second and fourth order error terms, and for a Wi 90° 
pulse this occurs at (/)2 = 2><pi and 0i = arccos(— 1/8), 
the well known BBl result 133] ■ The process can then 
be extended by using larger numbers of adjustable phase 
angles and seeking to remove higher order terms. 

With only two adjustable phases it is straightforward 
to tackle the problem analytically, seeking to set Taylor 
series coefficients to zero, but as the number of pulses is 
increased this method becomes intractable. Instead we 
have adopted a numerical approach, effectively replacing 



the solution of simultaneous non-linear equations with a 
minimisation problem by seeking to minimise the sum of 
squares of the first n coefficients: if this sum of squares can 
be reduced to a value indistinguishable from zero then the 
corresponding terms have been effectively eliminated from 
the series. (This approach is known to be unreliable in 
general }46i | , but seems to work fairly well in this particular 
case). Similar ideas have been explored elsewhere 30|,|40|. 
Our numerical searches show that adding a third ad- 
justable phase does not help; in particular the sixth order 
error cannot be eliminated. However adding a fourth ad- 
justable phase permits both the sixth order and the eight 
order terms to be eliminated with a sequence of eight tt 
pulses. This pattern continues throughout the range of 
searches we have conducted: the W„ sequence of 4n pulses 
with 2n adjustable phases 



7^</-l 71'02 



"^02,. '^<l>2 



7^02 ■^01 



(33) 



permits all the error terms up to order An to be removed. 
For the case of Wi there is only a single solution (the well 
known BBl composite pulse), neglecting the trivial vari- 
ant formed by negating all the pulse phases, but for W2 
there are two distinct solutions. These appear in theory to 
be effectively identical; in particular the tenth order coeffi- 
cient has the same size in the two cases. At W3 the search 
becomes computationally challenging, but once again two 
distinct solutions have been located. All these solutions for 
90° pulses are listed in Table [31 Preliminary searches indi- 
cate that at least one W4 sequence exists, but this has not 
been accurately located and it is not yet known whether 
any other sequences exist. 

This approach can of course be extended from 90° 
pulses to other angles; in particular it is possible to find 
phases for 180° pulses, as listed in Tabled These com- 
posite 180° pulses can then be converted into antisym- 
metric forms, and interestingly these appear to have the 
same property as Fi pulses, allowing them to be recur- 
sively nested. However this point is not explored further 
here. 



90° 


01 


02 


03 


04 


05 


06 


Wi 


97.2 


291.5 










W2 


84.3 


162.0 


345.5 


286.7 






W2 


132.3 


339.1 


26.4 


222.2 






W3 


22.0 


186.0 


89.3 


319.2 


178.2 


325.7 


W3 


79.9 


119.2 


257.3 


81.0 


308.4 


286.9 



Table 3: Composite 90° pulses: the Wn correction sequence contains 
4n pulses with 2n adjustable phases, and the phases (in degrees) 
listed permit all error terms up to order 4n to be eliminated when 
the correction sequence is combined with a naive 90° pulse. 



180° 



1 


104.5 


313.4 






2 


79.2 


193.4 


24.9 


307.5 


2 


130.0 


1.5 


56.7 


259.9 


3 


341.9 


147.9 


100.5 


355.9 



w,. 



207.3 339.4 
69.5 141.7 289.4 121.4 350.1 307.3 



Table 4; Composite 180° pulses: the phases (in degrees) listed per- 
mit all error terms up to order An to be eliminated when the Wn 
correction sequence is combined with a naive 180° pulse. 



8. Experiments 

Experiments were carried out in Oxford using a Var- 
ian Unity Inova system with a nominal ^H frequency of 
600 MHz to study a sample of HOD doped with GdCla 
to reduce the relaxation times. Experiments in Okayama 
used a homebuilt system with a ^^F frequency of 376 MHz 
to study a sample of perfluorobenzene dissolved in deuter- 
ated benzene. These two systems differ substantially in 
experimental detail, potentially allowing the effects of dif- 
ferent experimental errors to be distinguished. As our sys- 
tems should not suffer from significant off-resonance errors 
it is possible to introduce pulse rotation angle errors either 
by adjusting the pulse length (pulse-length errors) o r by 
adjusting the Bi field strength (pulse strength errors) |3l| . 
We chose to introduce deliberate errors by adjusting the 
pulse length, but as discussed below the system also con- 
tains underlying pulse strength errors. 

We begin by demonstrating the use of composite 90° 
pulses based on W„ correction sequences; for simplicity we 
simply show their performance as excitation pulses, but 
as Class A composite pulses the performance should be 
very similar for other initial states. The Bi field strength 
on the Oxford system was carefully adjusted to 25 kHz, 
corresponding to a nominal 90° pulse length of 10 fis, and 
then naive and composite pulses were applied using pulse 
lengths varying between 1 and 19 /iS. The signal intensity 
was determined by integration, and as the vertical scale is 
largely arbitrary all intensities were normalised by dividing 
them by the intensity from a naive pulse with e = 0. The 
results are shown in Fig. [51 with a vertical expansion in 
the lower panel; note that in this figure the lines simply 
join the experimental data points and are plotted to guide 
the eye. All experiments used the first of the the two 
choices for the W2 and W3 sequences; initial experiments 
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Figure 8: Experimental signal intensity (arbitrary units) as a func- 
tion of pulse length error e for naive (red circles), Wi (green trian- 
gles), W2 (blue squares), and W3 (black diamonds) 90° pulses used 
as excitation pulses. The lines simply join the experimental points 
and are plotted to guide the eye. The lower panel shows a vertical 
expansion of the upper panel. 



(data not shown) indicated that very similar results were 
obtained for the other two sequences. 

It is clear that the W„ composite pulses perform largely 
as expected, with wider error tolerance for higher values 
of n. However closer examination of the results, aided 
by the vertical expansion in the lower panel, shows two 
further effects. Firstly the overall signal intensity is larger 
with composite pulses than with naive pulses, and secondly 
all four plots lean to the right, showing higher intensities 
for positive values of e in comparison with the equivalent 
negative values. Both of these effects can be ascribed to 
additional pulse strength errors arising from Bi inhomo- 
geneity. The visible increase in overall intensity indicates 
that a significant fraction of the sample experiences a Bi 
field substantially different from the central value, while 
the lean to the right indicates that these regions experi- 
ence a smaller Bi field than the bulk of the sample. This 
behaviour is expected for traditional RF coil designs. 

Attempts were made to reduce these effects in exper- 
iments on the Okayama system by replacing the conven- 
tional sample tube by a small spherical sample in the cen- 
tre of the RF coil, which should reduce the RF inhomo- 
geneity over the sample. The results, shown in Fig. [9l 
indicate partial success. In particular the increase in over- 
all intensity is reduced, but a lean to the right remains 
visible. 

To make further progress we used narrowband n pulses 
to select signal from regions of the sample with a highly 
homogeneous Bi field strength. The double pulse field 
gradient spin echo (DPFGSE) sequence [47| is effective in 
selecting for spins which experience a 180° pulse at the 
centre of each spin echo. The classic use is to suppress 
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Figure 9: Experimental signal intensity (arbitrary units) as a func- 
tion of pulse length error e for naive (red circles), Wi (green trian- 
gles), W2 (blue squares), and W3 (black diamonds) 90° pulses used 
as excitation pulses with a small spherical sample. 



water signals by combining a hard tTx pulse with a soft 
■K-x pulse which selects water transitions, but it can be 
used more widely to select subsets of spins. Here we use 
an N2 composite tt pulse which only excites spins where 
the Bi field strength is close to its nominal value, and 
compare the results with those from an F2 pulse, which 
should pass a very wide range of field strengths. (This is 
simpler than comparing the N2 filtered data with spectra 
acquired without a selective filter, as the filtration process 
also causes signal loss due to the effects of diffusion |48J.) 

These filters were explored using a nutation sequence 
to measure RF inhomogeneity, as shown in Fig. [TUl The 
nominal Bi nutation frequency was adjusted to 25 kHz by 
adjusting the RF power until a 40 /.ts pulse corresponded to 
a notation of 360°, and the nutation time varied between 
and 127 /xs. The N2 data is well described by a sin- 
gle narrow Gaussian distribution of nutation frequencies, 
with a central value of 25 kHz and a width of 2.5 kHz, but 
the F2 data needs a second Gaussian component, centred 
around 20 kHz with a width of about 10 kHz. This low 
frequency broad component leads to the deviations from 
ideal behaviour seen in Fig. |S] and discussed above. 

Finally we show the effect of combining a W„ exci- 
tation pulse with a N2 filtration sequence. In this case 
the aggravating pulse strength errors arising from RF in- 
homogeneity should be well suppressed, leaving only the 
artificially imposed pulse length errors. The result, shown 
in Fig. [TT] confirms this, with the experimental data points 
lying very close to theoretical predictions. 

Note that in this plot the theoretical curves are sim- 
ply plotted directly on the same scale as the experimental 
data, rather than being optimised to fit this data, show- 
ing the very close match between theory and experiment. 
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Figure 10: Nutation experiments to measure RF inhomogeneity. In 
the upper panel the experimental signal intensity (arbitrary units) is 
plotted as a function of pulse length with F2 (red) and N2 (green) 
DPFGSE filtration sequences. The points show experimental values 
while lines show fitted curves. The N2 data is well described by 
a single narrow Gaussian component centred at 25 kHz, while the 
F2 data requires a two Gaussian fit, revealing a broad component 
centred around 20 kHz. The lower panel shows the distribution of 
nutation frequencies corresponding to these fits. 
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Figure 11: Experimental signal intensity (arbitrary units) as a func- 
tion of pulse length error e for naive (red circles), Wi (green trian- 
gles), W2 (blue squares), and W3 (black diamonds) 90° pulses used 
as excitation pulses followed by the use of a DPFGSE N2 sequence to 
select a small region with high Bi homogeneity. The points show ex- 
perimental values while lines show the theoretical curves. The lower 
panel shows a vertical expansion of the upper panel. 
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The remaining mismatch is most severe for the W3 exper- 
iment, and may partly reflect spin-spin relaxation during 
this long composite pulse. Remaining errors are likely to 
be due to imperfections in the phase-control of the RF 
source, and slow drifts in the power of the RF amplifier. 

9. Conclusions 

The ability to interconvert time symmetric and time 
antisymmetric tt pulses allows us to interpret previous re- 
sults with a common framework, and then to extend these 
results to produce new analytic families of composite tt 
pulses with unprecedented tolerance of pulse strength er- 
rors at little or no cost in sensitivity to off-resonance errors. 
This approach is currently confined to tt pulses, and can- 
not be applied to pulses with other rotation angles. It is, 
however, possible to use numeric methods to design sym- 
metric composite pulses for arbitrary angles, and several 
families have been located. Experimental implementations 
confirm that these pulses work very much as expected, but 
it is necessary to perform these experiments carefully in 
order to avoid confusion arising from multiple sources of 
pulse strength error. 
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